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ABSTRACT 

Based on long-slit infrared spectroscopic observations, it has been suggested that half of the 
carbon monoxide present in the atmosphere of comet C/1995 01 (Hale-Bopp) close to perihelion 
was released by a distributed source in the coma, whose nature (dust or gas) remains unidentified. 
We re-assess the origin of CO in Hale-Bopp's coma from millimeter interferometric data and a 
re-analysis of the IR lines. 

Simultaneous observations of the CO J(l-O) (115 GHz) and J (2-1) (230 GHz) lines were 
undertaken with the IRAM Plateau de Bure interferometer in single-dish and interferometric 
modes. The diversity of angular resolutions (from 1700 to 42 000 km diameter at the comet) 
is suitable to study the radial distribution of CO and detect the extended source observed 
in the infrared. We used excitation and radiative transfer models to simulate the single-dish 
and interferometric data. Various CO density distributions were considered, including 3D 
time-dependent hydro dynamical simulations which reproduce temporal variations caused by the 
presence of a CO rotating jet. The CO J(l-0) and J (2-1) observations can be consistently 
explained by a nuclear production of CO. Composite 50:50 nuclear /extended productions with 
characteristic scale lengths of CO parent L p > 1500 km are rejected. 

Based on similar radiation transfer calculations, we show that the CO v = 1-0 ro-vibrational 
lines observed in comet Hale-Bopp at heliocentric distances less than 1.5 AU are severely optically 
thick. The broad extent of the CO brightness distribution in the infrared is mainly due to optical 
depth effects entering in the emitted radiation. Additional factors can be found in the complex 
structure of the CO coma, and non-ideal slit positioning caused by the anisotropy of dust IR 
emission. 

We conclude that both CO millimeter and infrared lines do not provide compelling evidence 
for a distributed source of CO in Hale-Bopp's atmosphere. 

Keywords: Comet Hale-Bopp - Comets, composition - Radio observations - Infrared observations 
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Introduction 



About two dozen molecules have been detected in cometary atmospheres in addition to 
a wealth o f radicals, atoms and ions fo r med by photodissoci ation, photoionization and coma 



chemistry (jBockelee-Morvan et al 



2004 



Feldman et al. 



20041 ). It is generally agreed that most 



molecules found in the coma are released from the nucleus, being present in condensed form or 
as trapped gases in water ice. However, there are pieces of evidence that distributed sources 
(also called "extended sources") of gas are present in come tary atmospheres (see the reviews of 



Bockelee-Morvan and Crovisier 



2002 



Cottin and Fray 



20081 ) . Molecules released from the nucleus 



and from extended sources exhibit different radial distributions in the coma. Molecules produced 
from the thermal or photo-degradat ion of organic gr ains should display characteristic production 



curves versus heliocentric distance ( Fray et al. 



20061 ). Radial distributions not consistent with 



nucleus production were foun d for H2CO, CO, and OCS (e.g., 



Meier et al 



1993 



DiSanti et al 



2003 



Dello Russo et al 



1998). Strong heliocentric variations of the H2CO, HNC, and CS mixing 



(e. 



Biver et al. 


1999b: 


Frav et al. 


2006; 


Milam et al. 


2006; 


Lis et al. 


2008) 



20081 ). The production 



of radicals (CN, C2) from grains is also suggested. The study of extended sources of gases can 
provide important clue s to the chemical nature of refractory organic compounds in cometary nuclei 



rtCottin and Frav 



2008). For example, laboratory experiments show that the decomposition of 



polyoxymeth ylene, a polymer of formaldehyde, ca n be the source of monomeric H2CO in cometary 



atmospheres (jCottin et al 



2004; 



Frav et al 



20061 ). 



The first evidence for an extended CO production was obtained from measurements of the 



local CO density in comet IP/Hal 
neutral mass spectrometer (NMS). 



ey along the path of t 



Eberhardt et al 



re Giotto spacecraft acquired with the 



(jl987l ) concluded that only about one-third of 



the CO was released from the nucleus, the remainder being produced by an extended source in the 
coma. The total CO production rate in comet lP/Hall ev was estimated to be 11% that of water 



at the time of the Giotto flyby ([Eberhardt et a 



through a jet enriched in CO, was proposed by 



1987 ). An alternative view, that Giotto flew 



Greenberg and Lil (|1998l ) to interpret the NMS 
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observations. 



Constraining the radial distribution of CO in cometary atmospheres from telescopic 
observations requires significant spatial resolution and a productive comet. Long-slit spectroscopic 
CO observations near 4.7 /im performed at < 1.5 AU in comet C/1995 01 (Hale-Bopp) have 
suggested that 50 to 90 % of the CO in this comet was or iginating from a distributed source 



(Disanti et al 



1999; 



DiSanti et al 



2001 



Brooke et al 



2003). Similar observations performed in 



C/1996 B2 (Hyakut ake) indicate that CO was released almost entirely (~ 75%) from the nucleus 



(DiS anti et al. 



20031 ). CO was abundant in these two comets with a total CO production rate 



20% relative to water. 

We present here the study of the radial distribution of CO in comet Hale-Bopp from 
observations of the J (1-0) (115 GHz) an d J (2-1) (230 GHz) lines un dertaken with the IRAM 



Plateau de Bure interferometer (PdBI) ([Bockelee-Morvan et al 



2009). Observations were 



performed both in single-dish and interferometric modes with angular resolutions probing coma 
radii from 830 to 20800 km, which encompass th e cha racteristic sca l e (500 km) of the CO 



extended source identified by 



DiSanti et al. 



(2001 



) and 



Brooke et al 



(2003). Hence, they are 



appropriate for the study of the origin of CO gas in Hale-Bopp's coma. 

There is observ ation al evidence that th e infrared CO ro-vibrational lines observed by 



DiSanti et al 



(|200lh and 



Brooke et al 



(2003) in comet Hale-Bopp are optically thick. Optical 
depth effects make brightness radial profiles flatter than when optically thin conditions apply, 
and can lead to misinterpretations if they are not properly taken into account. We present a 
re-analysis of these lines using radiative transfer calculations. 

Both the radio and IR CO data show evid ence for the presence of a strong spiralling CO je t 



2009 



Brooke et al 



2003 1 



when comet Hale-Bopp was close to perihelion (jBockelee-Morvan et al 
This feature is considered in the present study. 

The observations of the CO radio lines are summarized in Sect. 12.11 the methods and models 
used for their analysis are described in Sect. 12.21 and the results are given in Sect. [3j The analysis 
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of the infrared lines is presented in Sect. 01 A summary follows in Sect. [5j 



2. Analysis of the CO millimetric data 



2.1. Observations : summary 



The detailed description of the observations of the CO J(l-O) (115 GHz) an d J(2-l 



GHz) lines undertaken i n com et Hale-Bopp with the IRAM PdBI can be found in 



Bockelee-Morvan et al 



Henrvl (|2003 ): 



(230 



(|2009l ). They were performed on 11 March 1997 from 4.4 to 13.8 h UT. 
Comet Hale-Bopp was at a heliocentric distance of 0.989 AU and a geocentric distance of 1.368 
AU. Both lines were observed simultaneously using different receivers and a spectral resolution of 
0.13 km s . At that time, the PdBI comprised five 15-m dishes. The interferometer was used in a 
compact configuration CI, with baselines (the spacing between two antennas), ranging from ~ 20 
to ~ 150 m. Most of the time, the comet was observed in cross-correlation (interferometric) mode. 
However, every one hour short-time (1 min on source) single-dish observations (also referred as to 
ON-OFF measurements) were performed. The observing mode was position-switching to cancel 
the sky background, with the reference position taken at 5 ' from the comet position. 

When all interferometric data are considered, the full width at half maximum (FWHM) 
of the synthesized interferometric beam is 2" x 1.38" at 230 GHz, and 3.58" x 2.57" at 115 
GHz. The FWHM of the primary beam of the antennas is 20.9" at 230 GHz and 41.8" at 
115 GHz. Given the geocentric distance of the comet, these angular resolutions correspond to 
projected distances of 1650 to 41600 km (diameter). Inter f erome tric maps of the CO J(l-0) and 



J (2-1) emissions are presented in 



Bockelee-Morvan et al 



(2009). The signal-to-noise ratio in 



the center of the maps ranges from 25 (J(l-0)) to 140 (J(2-l)). Line areas (i.e., flux densities 
integrated over velocity in units of Jy km s _1 , the lines being Doppler-resolved) measured on 
the overall single-dish (ON-OFF) spectra averaging all data (-Fsd)> and measured at the peak of 
interferometric maps (Fi n %), are given in Table [TJ For both ON-OFF and interferometric data, 
the uncertainties in flux calibration are at most 10% and 15% for the J(l-0) and J (2-1) lines, 
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respectively (jBockelee-Morvan et al. 



20091 ) . 



2.2. Methods of analysis 

The possible extended origin of CO can be investigated by comparing single-dish Fsd 
and interferometric Fi nt fluxes (Table [T]). Extended source production, if resolved out by the 
interferometric beam, would result in a ratio R = F^b /Fiat higher than the value expected for a 
pure nuclear production. On the other hand, opacity effects need also to be carefully investigated: 
being more important for smaller beams, they increase the -Fso/^int intensity ratio with respect 
to optically thin conditions. The measured intensity ratios i?i_o and R2-1 (the subscripts refer to 
the lines) are given in Table [TJ with error bars including uncertainties in flux calibration. 

The CO radial distribution can also be studied from the brightness distribution observed 
in the interferometric maps. We preferred to study the direct output of the interferometer, i.e., 
the complex visibilities measured by the antenna pairs which sample the Fourier transform of the 
brightness distribution in a number of points in the Fourier plane (also called uu-plane). This 
approach allows us to avoid uncertainties in map analysis related to: 1) the ellipticity of the 
interferometric beam which results from the anisotropic coverage of the Fourier plane; 2) flux 
losses expected for extended sources due to the non coverage of the w-plane at baselines shorter 
than 20 m. Visibilities (integrated over velocity) recorded by the ten baselines over the whole 
observing period were radially averaged in the uv-plane over intervals of uu-radius a of 5 m (the 
uv-radius corresponds to the baseline length projected onto the plane of the sky). The line area 
determined on the averaged single-dish spectrum provides the visibility at a = 0. 

The visibility from extended source emission is expected to drop faster with increasing a than 
for more compact sources. As shown in the Appendix, for a density distribution varying according 
to r~ 2 , where r is the distance to the nucleus, the amplitude of the visibility V is expected to vary 
according to u -1 , provided molecular excitation is not varying in the coma and the line is optically 
thin. For an optically thick line, the power index of the variation of V with a will be <— 1. For an 
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extended density distribution which follows a Haser daughter distribution, we computed that V 
oc <7 -2 , if the parent scale length L p is larger than the size of the interferometric beam and the 
line is optically thin. If molecular excitation is varying in the coma (e.g., as a result of variations 
of the gas kinetic temperature), V(<r) will deviate from these laws. The measured variation is 
oc ^-^±0.02 anc j a ^-o.ssio.oe^ f or ^ e J (2-1) and J (1-0) lines, resp ectively (slightly different 



values were derived from linear fitting by 



Bockelee-Morvan et al 



20091 ). 



These two approaches were used by 



Boissier et al 



([20071 ) to study the radial distributions of 



H 2 S, SO, and CS observed with the IRAM PdBI. 



2.3. Modelling 



Synthetic ON-OFF spectra and vis ibility spectra were c omputed using an excitation model 



Boissier et al 



([20071 ). We summarize here the main 



and a radiative transfer code, following 
aspects of the computations. 

The calculation of the population of the rotational levels of the CO molecules considers 
collisi ons with H2O and electrons, and IR radiative pumping of the v( l-O) CO vibrational 



band (Crovisier and Le Bourlot 



1983 



Crovisier 



1987 



Biver et al. 



1999b! ) . The model provides 



populations as a function of radial distance r, given a H2O density law with r. For simplicity, 
we assumed an isotro pic H2O coma. T 



a r = 2 x 10" 14 cm 2 



( Colom et al 



Biver et al 



re collisional CO-H2O cross-section is taken equal to 
1999 bj), and the H^O p r oduct ion rate is Qh 2 o = 10 31 s _1 



19991 ). As shown in 



Bockelee-Morvan et al. 



(|2009l ) . PdBI observations, including 



ON-OFF measurements, were mostly sensitive to CO molecules with rotational populations at 
local thermal equilibrium (LTE). The large size of the CO LTE region can be explained by the 
high collisional rates in this productive comet and the large radiative lifetimes of the CO J = 1, 2 
rotational levels. 

Two nominal gas kinetic temperature laws were investigated: 1) a constant gas kinetic 
temperature throughout the coma equal to 120 K, which agrees with temperature determinations 
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pertaining to the r ~ 10 000 km coma region (jBiver et al 



1999a; 



DiSanti et al. 



2001 



this 



temperature provides a best fit to the relative intensities of the CO J (2-1) and J (1-0) lines 



measured in single-dish mode at PdBI (jBockelee-Morvan et al. 



2009), assuming pure CO nuclear 



production and an isothermal coma; 2) a temperature law T var increasing with increasing 
cometoc entric distance, as expected from coma hydrodynamical models which include photolytic 



heating (jCombi et al 



1999). The long-slit infrared spectroscopic observations of CO and other 



molecules undertaken in March- April 1997 i n comet Hale-Bopp sugg e st typical gas kinet ic 



temperatures of ~ 80-90 K at r ~ 10 00 km (IMagee-Sauer et al 



adopted the temperature law used by 



Boissier et al 



1999; 



DiSanti et al. 



200l|). We 



(2007) to interpret the PdBI observations of 



sulfur-bearing species. It yields T < 70 K at r < 1000 km and T= 120 K at r > 10 000 km (see 
footnote in Table [2]). For a few calculations, we investigate an alternate temperature law T var ( a i t ) 



with T, 



var(alt) 



T var for r > 1500 km and T = 79 K for r < 1500 km. 



In most calculations, the CO density distribution was computed using simple steady state 
isotropic models. We considered outflow either at a constant velocity (i.e., Haser models) or at 
a velocity which increases with increasing cometocentric distance. Indeed, gas acceleration is 
expected in the r = 1000-10 00 km region prob ed by the PdBI in the Hale-Bopp coma, as a 



result of photolytic processes (jCombi et al 



on the gas expansion velocity v exp (e.g., 



1999). The shape of radio lines provides constraints 



Hu et al 



1991 



From t 



Biyeretal 



re wid th of radio lines of parent 



(|1999al ) derived a value equal to 



molecules (including CO) observed by single-dishes, 
1.05 km s -1 for the period pertaining to the PdBI CO observations. This value is representative of 
the gas velocity at cometocentric distances of typically the radius of the field of view and therefore 
pertains to radial distances ~ 10 000 km. The width of the CO lines in the interf erometric bea m 
sugge sts that v exp ~ 0.9 km s _1 at r ~ 1000 km. For comparison, the model of I Comb i et al 



( 1999l ) yields to an acceleration from 0.91 to 1.12 km s _1 in the r = 1000-10 000 km region, in 
good agreement with our velocity retrievals. For the models with constant velocity (v cons t), we 
adopted v exp = 1.05 km s _1 . For the models which consider the variation of the velocity (v vaT ), 
we assumed v ex p = 

0.9 + 0.15 x log(r[km]/1000) km s" 1 for r > 1000 km, 

^cxp — 0.9 km s 1 for 
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r < 1000 km. Figure Q] shows t 



he H?S 2sn-2i 



13 March 1997 with the PdBI (jBoissier et al 



line at 216.7 GHz observed in comet Hale-Bopp on 



20071 ). Synthetic line profiles computed with the v v 



velocity law are superimposed. The width and wings of both the ON-OFF and interferometric 
spectra are reproduced. The spectral asymmetries present in the central part of the observed H2S 
spectra are due to spatial asymmetries in the coma. CO and H2S are expected to expand at the 
same expansion velocity. We checked that the v var law consistently explains the width and wings 
of the CO J(l-0) and J (2-1) line profiles. 

In theory, cometary line shapes contain information on the radial distribution of the 
molecules. Daughter molecules are expected to exhibit double- horn line shapes when the size of 
the beam does not exceed much their production scale length. The missing flux in the center of the 
line depends on the ratio of these two quantities. However in comet Hale-Bopp, CO line profiles 
were found to be affected by the presence of a CO rotating jet which int roduced excess emission 



alternatively in the 



2009 



Boissier et al 



jlue and red part of the profiles during the rotation (jBockelee-Morvan et al 



20101 ) . Such complexity prevents detailed fitting of the CO line shapes. 



Extended distributions were investigated using the Haser formula for daughter species with 
the parent and daughter scale lengths L p and L^ and the gas velocity v exp as parameters. For 
simplicity, both L p and L^ (= Leo) were assumed to be constant throughout the coma, i.e., we did 
not consider the effect of the variation of the expansion velocity on the scale lengths. Therefore, in 
all calculations we as sumed Leo = ^exp//5cO! with u exp = 1.05 km s _1 and /3co = 7.5 x 10^'s -1 



( Huebner et al 



1992). Such approximation has a negligible effect on the results. Indeed, the CO 



scale length is much larger than the field of view of the single-dish measurements, so the CO 
column density within the field of view is not much affected by small inaccuracies in the effective 
CO scale length. The relative contribution of the extended and nuclear productions of CO is given 
by the ratio of the production rates Qqo'-Qqq, where the subscripts N and E refer to the nuclear 
and extended productions, respectively. The total production rate is Qco=Qco+^co- 

We also consider in Sect. 13.31 much complex distributions where a large fraction of the CO 
production is within a spiralling jet. As mentioned above, a strong CO jet was indeed identified 
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in the coma of comet Hale-Bopp from these interfero metric observations. The CO distribution is 



mode lled using tridimensional, time-dependent codes (jBockelee-Morvan et al 



201(1 ). 



2009; 



Boissier et al 



T he brightness d i stribu tion of the two CO lines was computed using the radiative transfer 



code of 



Biver 



Boissier et al. 



(|2007l ). which includes both self-absorption and stimulated emission (see 
19971 . for detailed explanations). ON-OFF synthetic spectra were obtained by convolving 
the brightness distributions with the primary beam pattern of the PdBI antennas, described by a 
2D Gaussian. Visibi lities were calculated from the Fourier transform of the brightness distribution 



( Boissier et al 



20071 ). Maps of the CO line emissions were synthesized, considering the coverage of 
the Fourier plane obtained during the observations. The peak flux densities in the modelled and 
observed maps, both obtained by inverse Fourier Transform, can then be directly compared. 



3. The CO radial distribution from millimetric data 

3.1. ON-OFF and interferometric flux densities 

Table [2] gives flux densities and intensity ratios RMod = -^sd / -Fint computed with Qco = 
2.1 xlO 30 s _1 and several Haser distributions describing the local CO density. Except for the 
first entry, corresponding to optically thin assumptions, all results were obtained with radiation 
transfer calculations. 

The value Qco = 2.1 xlO 30 s" 1 was chosen so that the J(l-0) and J (2-1) ON-OFF line 
intensities computed with T = 120 K, a constant velocity and a pure CO nuclear production match 
precisely the measurements. This value is consistent with the tot al CO production ra te (sum of 



2001 



The other 



native and extended distributions) inferred from IR spectroscopy (Di Santi et al. 
models considered in Table [2] yield ON-OFF line intensities consistent with the measured values if 
one takes flux uncertainties into account. 

To quantify the goodness of the models in explaining the observed intensity ratios i?obs (Ri-o 



-12- 



or i?2-i), we computed for both lines the quantity: 

X 2 = (#Obs " RModf/^Robs (1) 

where Ai?obs is the uncertainty on i?obs given in Table [T] and i?Mod is the intensity ratio given by 
the models. Results are plotted in Fig. [2j Models resulting in \ 2 > 1 i n either line are considered 
as non satisfactory. 



3.1.1. Models with pure nuclear production 

The J(2-l) line observations can be fully interpreted by a nuclear CO production, provided 
the opacity effects are taken into account. Opacity affects mainly the flux density of the J(2-l) 
line in the interferometric beam: e.g., a reduction by 37% is observed for the model with constant 
temperature and velocity, while ON-OFF J (2-1) flux and J(l-O) fluxes are almost unchanged 
(compare the first two entries in Table [2|). When line opacity is considered, the R2-1 values 
obtained with either of the temperature and velocity laws (~ 18-19) are consistent with the 
measured value of 17.8 ± 3.8. x 2 values close to zero are obtained (Fig. [2]). 

On the other hand, models considering a nuclear production and a constant gas kinetic 
temperature fail to explain the ratio i?i_o = 11.6dbl.8 measured for the J(l-O) line, yielding values 
20-35% higher than observed. The discrepancy obtained for Ri-o cannot be solved by invoking 
an extended CO production because this process affects the intensity ratio in the opposite way 
(Fig. [2]). However, a good agreement is obtained when the variable temperature law is used: the 
population of the J = 1 CO rotational level increases with decreasing T, hence producing an 
enhancement of the J(l-O) emission from the inner colder parts. The R2-1 ratio is less sensitive 
to the assumed gas temperature law because the increase of J(2-l) local emission in the inner 
regions is counterbalanced by more significant self-absorption effects for the emergent emission. 

The -Ri_o and R2-1 ratios are both sensitive to the assumed velocity law which affects the 
CO density radial profile. R±-o decreases when the expansion velocity in the inner coma decreases, 
the outflow velocity at r > 10 000 km being kept constant. R2-1 displays a non monotonous trend 
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with decreasing inner coma velocity v exp (inner), slowly decreasing for 0.9 km s 1 < v exp (inner) 
< 1.05 km s" 1 (Table E}, and then increasing for v exp (inner) < 0.9 km s . We found that it is 
not possible to explain both Ri-o and R2-1 measured ratios with a variable velocity law and a 
constant gas kinetic temperature. 

In summary, the measured intensity ratios R\-o an d R2-1 can be explained with a CO nuclear 
origin when the variation of the gas kinetic temperature with distance to nucleus is considered. 



3.1.2. Models with extended CO production 



Let us first compare th e measurements to model results obtained with pure CO daughter 
Haser distributions. Indeed, 



Brooke et al 



(2003) interpreted long-slit infrared observations of CO 



in comet Hale-Bopp, with ~ 90% of the CO being released from a source with L p = 5000 ± 3000 
km. We investigated CO parent scale lengths up to L p = 5000 km. Such CO distributions with 
L p > 500 km do not fit the data (Fig. [2D- In other words, both the CO J(l-0) and J(2-l) data 
exclude that CO is produced solely by an extended source, unless this source releases most of the 
CO gas in the near-nucleus coma unresolved by the interferometric beam. Results obtained with 
Qqq'-Qqq in ratio 10:90 and L v = 5000 k m do not rep r oduce the measurements (Table [2]), and 



thereby do not support the conclusions of 



Brooke et al. 



|2003j). 



DiSanti et al. 



(|200ll ) interpreted the radial profiles of CO infrared emission by the contribution 
of a distributed source to 50 % of the total CO production. Therefore, we also investigated 
composite nuclear:extended productions with CO and CO parent production rates in ratio 50:50 
(Table [21 and left part of the plots shown in Fig. [2]). Satisfactory fits to both lines can only 
be obtained for models with the variable temperature law: the i?i_o and R2-1 ratios can be 
reproduced within 1-a providing L p < 1000 km (model with constant velocity) or L p < 2000 km 
(model iVar)- One can consider that the condition X1-0 + X2-1 < 1 should be satisfied as the 
J(l-0) and J(2-l) measurements are independent. Then, models with L p < 1500 km (v vauT law) 
and L p < 800 km (v cons t law) provide a satisfactory fit to the PdBI CO data based on the ratios 
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R of single-dish to interferometric fluxes. Since the adopted v var law is consistent both with line 
shape measurements and gas dynamics calculations (Sect i2^3|) . we conclude that composite 50:50 
nuclear /extended productions with L p > 1500 km are excluded by the CO millimeter data. 

It is worth noting that, using the models with, the tVar and T var laws, and extended 
productions fitting the data, we infer total CO production rates within 5% of the nominal value of 
2.1 xlO 30 a -1 . 

Still, we assumed isotropic CO density distributions. In Sect. 13.31 we discuss how the 
Fsr>/Fi ni ratio is affected by the presence of jets. 

3.2. Radial evolution of the visibilities 

In Figs [SHU are plotted the visibilities V as a function of uu-radius a expected for different 
models, together with the measurements. The power index of the variation of V(cr) with a 
measured on the modelled curves is given in Table [2j The least-squares fit to the computed 
visibilities was performed on an extracted set of visibilities representative of the a sampling of 
the observations. Weighting factors corresponding to the uncertainties in the measured visibilities 
were considered when fitting the extracted set of computed visibilities. This method is appropriate 
for best comparison with the data, since the model results show that the visibility curves do not 
follow a simple power law (Fig. . 

From visual comparison between the modelled and observed visibility curves (Figs [SHU, we 
see that : 

1. Both lines exclude pure CO production from an extended source, unless the characteristic 
parent scale length of the distributed source is well below 500 km; 

2. Composite nuclear/extended productions with L p = 5000 km or 2000 km do not explain the 
observations. 
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3. The best fits are obtained for a nuclear production. 

From a least square fitting of the visibility curves, including the ON-OFF point, we found 
that composite 50:50 nuclear /extended productions with L p > 1500 km are not statistically 
favored when the v wai law is considered. In other words, we obtained the same results as derived 
from the intensity ratios. 

Let us now compare the power indexes of modelled visibility curves to the measured values. 
For the J(l-O) line, there is a good agreement (within 1.3-sigma) between the measured value and 
that given by the model with a nuclear CO production, and T var and v var laws. On the other 
hand, the same model predicts a power index of —1.40 for J(2-l) whereas the measured value 
is —1.27 ± 0.02. The steep decrease of the J(2-l) visibility curve obtained with the T var law at 
large un-radii (> 70 m, see Fig. U]) is caused by large self- absorption effects in the inner cold 
coma observed by the interferometric beam (r < 850 km) where T < 65 K. At 230 GHz, uu-radii 
between 70-150 m probe cometocentric distances (oc A/cr) to which 115 GHz measurements are not 
sensitive. Unlike J(2-l) visibilities, J(l-0) visibilities are not influenced by the coma temperature 
at r < 1500 km. Figure [5] shows that the J(2-l) data at large ttf-radii can be better fitted using 
the alternate temperature law T mr ( a ^) introduced in Sect. 12.31 for which T = 79 K at r < 1500 
km and T = T var at r > 1500 km. The power index of the visibility curves computed with this 
alternate temperature law are —1.35 and —0.98 for J(2-l) and J(2-l) (Table [2]), respectively, in 
reasonable agreement with the measurements. We plot in Fig. [6] curves of \ 2 values obtained with 
the T var ( a itj law, which are analogous to those shown in Fig. [2j 

The model with composite 50:50 nuclear /extended production and L p = 5000 km closely 
reproduces the .7(2-1) measured power index, but not the flux ratio R2-1 (Table [2]). It is unable 
to explain either the power index of the J(l-0) visibility curve and R\~q- This shows that best 
constraints on the radial distribution of CO are obtained from the comparison of ON-OFF and 
interferometric flux measurements. 
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3.3. Models with anisotropic CO production 



A question that may arise is the extent to which our approach is appropriate in presence of 
jets in the coma. A spi ralling CO jet comprising ~ 3 % of the total CO production was identified 



from these PdBI data (Bockelee-Morvan et al. 



20091 ). This jet, found to be issued from a low 



latitude nucleus region and to be ~ 40° wide, caused temporal modulations in the visibilities 
recorded by each pair of antennas, which are still visible in the radial and time averages shown in 
Fig. E] and H 

The question has been thoroughly studied for spiralling jets of various streng ths and latitu des 



in the assum ption of negligible optical depth and a constant coma temperature (jHenrv 



200J) 



Henrvl (|2003l ) used a simple description of the Hale-Bopp CO coma represented by the combination 
of an isotropic coma and a spiralling conical jet. These 3-D time dependent simulations show 
that, for nearly equatorial jets, the power index of the visibility curve remains close to —1, except 
for collimated jets where the slope can reach —0.8 to —0.7 depending on the jet strength. An 
opposite trend is observed for high-latitude je ts (see results obtained for optically thin conditions 



in Figs 20, 21 of 



Bockelee-Morvan et al 



2003 ) . These trends are explained by the sampling time 
of the measurements (2/3 of the nucleus rotation, i.e., almost one full rotation), the orientation 
of the different baselines with respect to the jet direction, and the shape of the modulations. The 
shortest baselines probed more particularly gas emissions along the rotation axis , whil e the longest 



baselines were sensitive to equatorial jets (see Figs 18 of 



Bockelee-Morvan et al 



20091 ). Hence, in 



presence of collimated equatorial jets, the slope of the visib ility curve decreases. Simulations we re 



2003; 



Bockelee-Morvan et al 



not p erformed for the J(l-0) line in presence of a CO jet (jHenrvi 
2009). However, we do not expect conclusions different from those obtained for J(2-l) since 
J(l-0) and J(2-l) were observed simultaneously, i.e., with the same baseline orientations. Less 
contrasted visibility modulations and slope variations are even expected for J(l-0) in presence of 
jets, since the spatial resolution is degraded at 115 GHz. 



In the simulations of the CO PdBI data performed by 



Bockelee-Morvan et al 



((2009J), optical 



depth effects are considered and the kinetic temperature is assumed to be equal to 120 K. When 
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the jet parameters that provide the best fit to the data are considered (jet model (3)), the 
power index of the J (2-1) vis ibility curve agrees with that measured on the observed visibilities 



(Bockelee-Morvan et a! 



2009|). 



We also investigated to which extent the flux ratio R21 is affected by the presence of the CO 
jet. A value of 18 .4, consistent with obse r vation s, is derived from the synthesized interferometric 



maps obtained by 



Bockelee-Morvan et al 



(2009) with their jet model (3). Other m odel simulatio ns 



are worth to consider. In contrast to the simple descrip tion of the coma made by 



and 



Bockelee-Morvan et al 



|2009j) 



Boissier et al 



Henrvl (|2003h 



(|2010i ) simulated the CO PdBI observations of 



comet Hale-Bopp using a tridimensional time-dependent gas-dynamical model of the coma. An 
increased CO production from a localized region on the nucleus surface is considered to produce 
the observed CO spiralling structure. R21 values in the range 18-18.6 are inferred from their 
solutions 3-5 which provide the most satisfactory fits to the data. These values are comparable to 
the values found under the assumption of isotropic CO production (Table [2]) and to the observed 
value. 



4. The CO radial distribution from infrared data 

4.1. General considerations 

The millimetric and infrared CO data lead apparently to inconsistent results concerning 
the origin of CO in Hale-Bopp's coma. Based on radiative transfer calculations of CO v = 1-0 

1 — 1 1 — n 

ro- vibrationa l line e mission at 4.7 fim, we argue that the analyses of iDiSanti et all (|200ll ) and 



Brooke et al 



(2003) were hampered by an improper account of the opacity of these lines. 



Qualitatively, opacity effects are an attractive explanation: 

» Being more important for the central pixels aimed at the inner coma, they make line 
brightness radial profiles to be flatter than when optically thin conditions apply, and to 
mimic profiles for species with extended production; 
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They are all the more important sinc e the CO production rate is large and the geocentric 



distance is low. 



DiSanti et al. 



(2001 



J) found that in comet Hale-Bopp the CO extended 
source was present for < 1.5 AU (Qco/A in the range 5-16 x 10 29 s _1 AU _1 ) and 
inactive for > 2 AU (Qco/A < 1 x 10 29 s _1 AU _1 ). Interestingly, comet Hyakutake was 
also a CO-rich comet and the CO infrared data show ed extended brightness profiles when 



Qco/A was 1 to 5 10 29 s^AlT 1 (jDiSanti et al 



20031 ). Since the optical thickness increases 



with decreasing gas kinetic temperature, one expects opacity effects to be relatively more 
prominent in the CO brightness profiles of comet Hyakutake where the gas temperature 
was measured to be smaller. We are not aware of any other comet observed by means of 
high-resolution infrared spe ctroscopy showing evidence of an extended CO distribution 



(DiS anti and Mumma 



20091 ). So far, for all other comets studied to date in the infrared, 
Qco/A was much lower than 10 29 s AU . 

• Comparison of the intensity of lines emitted by a common ro-vibrational level provides 
direct evidence for significant opacity effects in the lines of sight passing close to the nucleus. 
Under optically thin conditions, the ratios of P2 to R0 (upper rotational level J u = 1), and 
of P3 to Rl (J u = 2) line intensities are equal to 2 and 3/2, respectively (i.e. (J u +1)/J u ). 
Smaller values are observed, consistent with the higher absorpti on coefficients of the P2 

). The 



and P3 lines as compared with those of R0 and Rl, respectively ( DiSanti et al 



2001 



radial brightness profile of R0 is less extended than those of R l, R2 and R3, whic h can be 



explained by the lower R0 absorption coefficient (see Fig. 13 of lBrooke et al. 



20031 ). 



Opacity effects affect the fluorescence excitation of the CO v = 1 vibrational band in the 
inner coma, the solar photons being absorbed along their path. Opacity effects also result in 
self-absorption of infrared photons emitted by nearby CO molecules (radiation trapping), which 
counterbalances the reduced direct excitation by solar photons . Last, but no t least , the received 



DiSanti et al. 



(|200ll ) considered 



radiation can be severely decreased by self-absorption effects, 
the attenuation of the solar pump and computed effective g-factors accounting for optical depth 
effects. They did not studied quantitatively the fate of emitted photons, and assumed optically 
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thin lines. Using effective g- factors, their conclu sions remained essen tially unchanged concerning 
the extended source of CO in comet Hale-Bopp. 



DiSantietal 



(|2003l ) found that the attenuation 
es observed in comet 



Brooke et al. 



(120031 ) 



of the solar pump explains in large part the CO extended brightness profi 
Hyakutake, especially those obtained end of March 1996 at A = 0.1 AU. 
considered optical effects both in the solar pump and in the emitted radiation to explain data 
obtained on 5 March 1997, and found spatial profiles only slightly wider than in the optically thin 
case. 

We note that, despite CO being less abundant than water, its v = 1-0 ro-vibrational lines are 
expected to be much thicker than the H2O lines from hot bands observed in the infrared. Photons 
emitted from a ro-vibrational transition (v\ i)— >(v" , j) can only be re-absorbed by the (V', j) level 
to excite (y\ i). In cometary atmospheres, most of the molecules are in their ground fundamental 
vibrational state. Hence, self-absorption is much less important for H2O hot bands, which connect 
weakly populated high vibrational states, than for the C O v = 1-0 band. Using the infrared 



excitation model of 



Bockelee-Morvan and Crovisiei ( 



H2O lines observed in comet Hale-Bopp (jDello Russo et al. 



1989), we computed that the op acity of ^3-^2 



2000; 



Brooke et al. 



2003 ) is typically 



5-6 orders of magnitude lower than the opacity of CO v = 1-0 lines at 1 AU from the Sun. For 
example, we can compare the P3 (J = 2 — > 3) line of CO v = 1-0 (2131.6 cm" 1 ) to the 2i2 — » 
3o3 line of H2O 1^3— 1/2 (2003 cm -1 ). The line absorption coefficient, and therefore the opacity, 
is proportional to the Einstein coefficient for absorption Bji (oc Aij, the Einstein coefficient for 
spontaneous emission) times the population of the lower level of the transition rij. Assuming 
that the rotational levels of the ground vibrational state are populated according to a Boltzmann 
distribution at 80 K, for the CO line we have = 20.3 s _1 and nj = 0.14, while for the H2O line 
= 1.14 s" 1 and nj = 7.0 10~ 6 . 
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4.2. Modelling 



We performed detailed calculations of CO v = 1-0 line intensities solving exactly the 
radiative transfer by integration through the coma. The source functions at each radial distance 



19871 ). This method 



were determined using the escape probability formalism (|Bockelee-Morvanl 
allows one to solve the statistical equilibrium equations for the level populations decoupled 
from any radiative transfer. It accounts for absorption of nearby emitted photons, stimulated 
emission, and for the attenuation of the solar flux by a correction factor to the excitation rates. 



Zakharov et al 



(2007) showed that the escape probability formalism is appropriate to account 



for radiation trapping in cometary atmos pheres. Escape probab ilities were computed assuming 



sp herical expans i on at constant velocity (jBockelee-Morvan 



by 



Brooke et al. 



19871 ). The same approach was used 



(2003). To compute fluorescence excitation of the v = 1 vibrational band, the 
population distribution in the ground vibrational state is needed. We used the excitation model 
for rotational emission described in Sect J2.2i For completeness, collisional quenching of v = 1-0 
fluorescence was considered assuming a collisional cross-section for de-excitation equal to 2.5 x 
10~ 16 cm 2 , H2O being the collision partner. However, this process was found unimportant for 
comet Hale-Bopp close to perihelion, except within a hundred of kilometers from the nucleus 
surface. Because of Doppler shifts, the optical depth experienced by the CO emitted photons 
depends on the velocity distribution and radial profile. As for the radiative transfer calculations of 
CO rotational line emission, constant or variable velocity radial profiles were assumed. The local 
velocity dispersion was described by thermal broadening. 

Two radiative transfer codes were developed. In the first one, which results are presented 
in Sect 14.41 isotropic Haser distributions were considered. The second one is an adaptation of 



the code of 



Bockelee-Morvan et al. 



(|2009l ) used to interpret the CO interferometric radio data 
(see Sect J3.3p . The CO coma is modelled as the sum of an isotropic contribution and a spiralling 
jet. The gas velocity and tempe rature are assumed to be the same in the two components (see 



discussion in 
Sect. 1431 



Boissier et al. 



2010|). Calculations performed with the second code are presented in 
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To study the extent to which the spatia 
expected for direct release from the nucleus, 



4.3. Q-curve analysis 

distribution o f com pounds deviates from that 



DiSanti et al 



(|200ll ) developed the method of 



Q-curves. An apparent production rate Q app is derived from the intensity Fi ine (x) measured at a 
specific location x along the slit, under the assumption of spherical outflow from the nucleus at 
constant velocity : 

hcug Une f(x) 

where gu ne is the line fluorescence g-factor, and /3 is the CO photodissociation rate. The quantity 
f(x) represents the fraction of the total number of molecules in the coma within the region 
sampled at offset x. Q-curves trace the evolution of Q a pp( x ) a s a function of x (given in arc sec or 
in km with respect to the nucleus position). For direct release from the nucleus and optically thin 
conditions, Q-curves reach quickly a terminal value equal to the true molecular production rate. 
Qapp values measured at small offsets from the nucleus underestimate the true pr oduction rate 



as th e measured intensities are here affected by slit losses due primarily to seeing (jDiSanti et al 



2001 



). The ratio between the terminal and on- nucleus Q a pp values is larger for extended brightness 
distributions than for more compact distributions. An extended distribution can be identified by 
comparing with Q-curves obtained for dust or species of purely nuclear origin observed in similar 
seeing conditions. To examine how optical depth effects may have affected the CO Q-curves, it is 
then important to take into account slit losses in our modelling. 



4.4. Results with isotropic distributions 

In this section, we present results obtained under the assumption of an isotropic CO coma. 
Figure [7] shows Q-curves computed for conditions corresponding to the data acquired on December 
11, 1996 (r h = 2.02 AU, A= 2.83 AU, Qco = 3.0 x 10 29 s^ 1 ) and January 21, 1997 (r h = 1.49 
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AU, A = 2.2 AU, Qco = 1-1 x 10 s )• Carbon monoxide is assumed to be solely released 
from the nucl eus. The CO r o -vibra tional lines used for generatin g these Q-cury e s are the same as 



those used by iDiSanti et al 



(120011 ) in their Fig. 7. According to IDiSanti et al 



POOJ ), the data 



of December 1996 are consistent with direct release of CO from the nucleus, whereas those of 
January provide evidence for production from an extended source. We convolved the modelled 
brightness distribution of CO v = 1-0 emission lines by a 2-D Gaussian representing the point 
spread function (PSF), and computed expected fluxes in 1 x 1 " boxes along the rad i al dir ection, 



for direct comparison with the data of comet Hale-Bopp presented by IDiSanti et al.l (|200ll ) . We 
assumed a seeing equal to 3" and 2.1" for modelling the December and January data, respectively. 
These values provide reference dust Q-curves (see below) consistent with those observed. The 
gas kinetic temperature was taken equal to 60 K (Dec. data) and 80 K ( Jan. data), based o n 



the temperature measurements in the outer coma from the CO IR lines (Di Santi et al. 



2001 



)• 



Calculations shown in Fig. [7] were performed under the assumption that the gas expansion velocity 
v exp is constant through out the coma and equal to 0.9 km s _1 (Dec. data) and 1.0 km s -1 (Jan. 



data) ( Biver et al 



2002 ). 



For each date, two Q-curves were generated: 1) a reference Q-curve (empty circles) where the 
lines were assumed to be optically thin, collision quenching of fluorescence emission was neglected, 
and the rotational level populations in the fundamental were assumed to be at LTE; this reference 
Q-curve is also representative of the expected Q-curve from dust, providing dust emission falls off 
as 1/p, where p is the projected distance from the nucleus; 2) the second Q-curve (filled circles) 
was generated from the intensities computed with the full model. Apparent production rates were 
derived using Eq. [2] and g-factors at T = 60 K (Dec. data) and 80 K (Jan. data). The comparison 
of these Q-curves shows that the CO ro-vibrational line emission was severely affected by optical 
depth effects in January 1997 (Fig. [7J. The terminal Q app values for the scaled refe rence Q-curve 



and C O Q-curve are in ratio Rq = 0.56. This is close to the measured value of 0.49 (jDiSanti et al 



200ll ). showing that the opacity of the lines explain in large part the extended nature of the spatial 
brightness profiles. At the largest offsets plotted in the figure, Q app is slightly below the true CO 
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production rate because opacity effects are still significant. As for the December conditions, the 
model predicts that opacity effects should be discernible on the Q-curve: the scaled reference 
Q-curve and CO Q-curve are predicted to be in the ratio Rq = 0.83. Figure [8] plots the mean line 
opacity of several R and P lines for the conditions of early March 1997 (A = 1.47 AU, Qco = 2 
x 10 30 s -1 , T = 90 K). The mean line opacity is a weighted-average of the opacities in spectral 
channels covering the Doppler-resolved line, where the weight is the intensity in the spectral 
channels. The opacity of the R6 line (4.9 for early March) drops to 2.2 and 0.5 for mid-January 
and mid-December, respectively. Note that the IR CO lines are generally much thicker than the 
radio lines (Fig. [5]). 



As discussed by 



DiSanti et al 



(|200ll ). a reduction of the opacity of the CO infrared lines may 



be expected if one considers gas acceleration. To investigat e this effect, we solv ed the radiation 



transfer with gas accelerations consistent with the model of 



Combi et al 



((19991). In the law used 

to interpret the December data, the gas velocity increases from 0.81 to 0.92 km s _1 for radial 
distances between 1000 and 20 000 km. As for January, it increases from 0.81 (1000 km) to 1.03 
km s -1 (20 000 km). The results are shown in Fig. [9j A small reduction of the opacity effects is 
obtained. Rq values are increased from 0.83 to ~0.90 and from 0.56 to 0.60, for December and 
January, respectively. We consider that the December data are satisfactorily explained by the 
model. 

A 20% discrepancy between the computed and measured Rq values remains for the January 
data, which indicates that the observed spatial profiles are slightly more extended than expected 
for a pure C O nuclear origin. Cal culations of brightness profiles applying to the February to May 



1997 data of 



DiSanti et al. 



(|200ll ) show that a discrepancy of 20-30% is systematically present. 
Figure [TU] shows Q-curves for January conditions generated with combined 50:50 nuclear /extended 
distributions and CO parent scale lengths L p from 1500 to 5000 km. Rq decreases from 0.56 to 
0.45 for L p increasing from km (i.e., nuclear production) to 5000 km. Considering the reduction 
of opacity effects caused by gas acceleration, the best fit to the January data is obtained for 



L p = 5000 km. 
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While we cannot exclude a contribution of a distributed source of CO, its presence is not 
so compelling. Indeed, spatial profiles of water obtained with the same high resolution infrared 



instrumentation were also found to be broader t 



lan dust profiles and pr ofiles computed for a p 1 



brightness distribution convolved with the PSF (jDello Russo et al. 



2000) • The differences between 



water and dust /PSF Q-curves are comparable to or larger than the differen ces between observed 



Dello Russo et al 



(|2000i ) discussed several 



and computed (pure nuclear origin) CO profiles, 
possible explanations for the different behaviors of water, dust and PSF Q-curves and could not 
conclude for the contribution of a distributed source of water. In the next section, we argue that 
the complexity of the CO and dust Hale-Bopp comas should be considered to fully explain the CO 
spatial profiles. 



4.5. Model with anisotropic CO distribution 



As mentioned earlier, close to perihelion the coma of comet Hale-Bopp was structured by 
a strong CO spiralling jet or iginating from a low latitude region in the northern hemisphere 



( Bockelee-Morvan et al 



2009). In contr ast, the dust coma was dominated by a high latitude 



northern jet producing repetitive shells (jJorda et al 



1999 



Vasundhara and Chakrabortv 



i99a) 



This northern dust jet introduced a significant departure between the nucleus position and the 
center of brightness in visible images. Astrometric positions obtaine d around perihelion were 



all off by ~ -3" in Declination with respect to the nucleus position (jBoissier et al.l 



20071 ). The 



observations of the CO IR lines were performed with t he slit centred on the peak signa" 



continuum and aligned along the East- West direction (jDiSanti et al 



2001 



Brooke et al 



in the 



200J) 



Hence, they were likely not centred on the nucleus position. 

These spatial East- West profiles presented complex asymmetric shapes, with the maximum 
of the em ission frequently peaking at 1-2 x lO 3 k m (typically 1—1.5") West or East of peak dust 



emission (DiSanti et al. 



2001 



Brooke et al 



2003). These complex shapes might be related to the 



CO rotating jet. 
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We performed 3D radiative trans fer calculations for CO infrared lines using the coma 
model of iBockelee-Morvan et al.1 (|2009l ). We present results obtained with their jet model (3) 
for the geometric conditions of early March 1997, Qco = 2 x 10 30 s^ 1 and T = 90 K. Several 
sets of I R CO data have been acquired during this period (namely on March 1.08, 3.80 and 



5.99 UT DiSanti et al. 



2001 



Brooke et al. 



20031 ). which is close to the date (11 March) of the 



interferometric CO observations used to constrain this jet model. The jet is defined by its 
half-aperture (vl/ = 18.3°), latitude {I = 20°), and the fraction of CO molecules released within 
the jet (f co = 35.5%). In early March 1997, the spin axis had an aspect angle U ~ 77 ° anc . 



Jordaetal 



a posi tion angle on the sky pa^ ~ 205° based on the spin orientation derived by 
(|1999l ). Figure [TT1 shows the expected brightness distribution of the R6 CO line as a function of 
time over one nucleus rotation period. The jet structure is aligned along the RA axis towards 
West during about half of the period. The brightness distribution is m ore extended along the 



j et. U sing the longitude of the jet at the nucleus surface determined by 



Bockelee-Morvan et al 



2009) for 11 March UT a nd a synodic rotation period of 11.31±0.01 h (jFarnham et al 



Jorda and Gutierrez 



1999; 



2000), we inferred that the brightness distributions for 1.87, 3.80 and 5.99 



March UT should be approximatively those shown in panels 3, 4, and 11 of Fig. [TT1 respectively. 
Next, we will also consider panel 10 for comparisons with the 5.99 March UT data. 

The 1.87 March UT spatial profile peaks East of the continuum pe ak and falls off less rapidly 



with increasing p towards West than towards East (jDiSanti et al 



2001 



). This is consistent with 



spatial profiles extracted from panel 3 (Fig. [T2]) . As for the spatial profiles of 3.80 and 5.99 Ma rch 



Brooke et al. 



(|2003h 



UT, the sign of the offset would match the data if the directions given in 
were reversed. On 3.80 March the peak intensity is observed westward the continuum peak and 
is more extended towards East, while the model predicts the opposite (Fig. [T2]) . This is puzzling 
because the time between the March 1.87 and 3.80 UT observations is 4.1 periods, so we expect 
the profiles to exhibit the same shape. For 5.99 March UT, the peak intensity is towards the East, 



but predicted towards West (panel 11) . T. Brooke 



possible that the orientations given in 



Brooke et al. 



person al communication) comments that it is 



(2003) are in error though he does not have 
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firm evidence. 

Figure [T2] displays extractions parallel to the East- West direction for several negative offsets 
JDec in Declination. The peak intensity decreases when the offset increases, so that the spatial 
profile broadens. The position of the centroid of the emission reaches 1 to 3.5" for the times 
corresponding to panels 3, 4 and 10. Very broad profiles are obtained for panel 11 (expected to 
represent 5.99 March data), because the jet is nearly in the plane of the sky and aligned along 
East- West at that time (Fig. [TT]) . 

The intensity at 12" normalized to the peak intensity is in the range 0.12-0.20 at <5Dec = ", 
0.13-0.25 at <5Dec=-1.2", and 0.15-0.30 at <5Dec=-2.4". Values as high as 0.3 (<5Dec=-l") and 0.5 
(<5Dec=-2.5") are inferred in the West side of extractions from panel 2. For comparison, on the 
spatial profiles recorded in March- April 1997, this quantity generally ranges between 0.15-0.30, 
with extreme values of 0.10 and 0.4 being also observed. Therefore, the extended appearance of 
the CO brightness distribution could also be partly related to the positioning of the slit. Detailed 
modelling of the complex shapes of the CO IR radial profiles would be challenging. They are 
strongly sensitive to the .jet char a cteris tics and nucleus rotation properties. The time-dependent 



model of 



Bockelee-Morvan et al 



(|2009l ) provides a very simplistic description of Hale-Bopp's 
rotating CO coma, and does not explain in detail the interferometric CO maps. More sophisticated 
models would be required, and this is beyond the scope of this paper. 

So far, we did not discuss how the absolute intensities given by the model compare to 



Brooke et al 



(120031 ') . as 



the measurements. This exercise is possible for the data published by 
intensities are listed in the paper. Table [3] lists measured R2 and R6 line fluxes. Fluxes in a 
2"xl" inner box centred on the brightness peak (Fi nner ), and East- West averages of 5"xl" boxes 
starting 3" from the flux peak (F ou t er ) are given. Results from model calculations are also given 
in Table [3] for three slit offsets in Declination (<5Dec = 0, -1.2, and -2.4"). To compare to the 
5.99 March data, we show results obtained using the radial profiles from panel 10 instead of 
those from panel 11, as the latter are much broader than those observed. The measured line 
flux ratios F ou t er / Fi nner (~ 1.04-1.09, Table [3]) can be explained providing the slit was offset by 
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1-2" in Dec from the nucleus position. The intensities measured on 3 March suggest that the CO 
production rate wa s higher than 2 x 10 30 s -1 at that time. Using simple scaling, we infer a value 



consistent with thelBrooke et al. 



(|2003l ) determination of 2.6 x 10 30 s 1 based on the line fluxes 



at large projected distances (outer box). The intensities measured on 5 M arch suggest a total CO 



production rate of ~ 1.4 x 10 30 s , consistent with the value derived by 



Brooke et al 



(2003). 



Note that the IR H2O observations perform ed on 5 March also suggest that comet Hale-Bopp was 



less active on that day (jBrooke et al 



2003) 



Table provides also the intensities of the R2 and R6 lines obtained using a steady state 
isotropic model for the CO density distribution. The flux ratios F ou t er / Finner are, for most entries, 
slightly lower than the values obtained with the anisotropic model (i.e., the model with the jet 
yields East- West averaged radial profiles only slightly broader than for an isotropic coma). We 
computed that, for a seeing of 1.7, the expected F outer / Fi nner ratio for isotropic and optically 
thin conditions is 0.3 3, consistent with v alues measured for CH4 and H2O in early March for 



similar seeing values ([Brooke et al 



2003). Hence, the broad spatial distribution of the CO IR line 
emission results mainly from optical depth effects and, at a lesser extent, from the anisotropy of 
the CO coma and non-ideal slit positioning. 



Our conclusion is inconsistent with t he one obtained by 



modelling, but isotropic CO distributions. 



Brooke et al 



3rooke et al 



(2003) using similar 



(|2003l ) analysed the radial profile of the 



R6 line observed on 5 March, and were able to fit the flux in the inner box Fi nner with Q{CO) = 4 
x 10 29 s _1 and T = 100 K. With such a low CO production rate, spatial profiles were only slightly 
broader than in the optically thin case. Using the same parameters, our model yields a value for 
Finner which is 40% lower than the observed value, and a F outer / Fi nne r ratio equal to 0.43, i.e., 
barely higher than the optically thin value. Our model requires a higher CO production rate to 
fit the observations, yielding more pronounced opacity effects. We were not able to elucidate the 
origin of the discrepancy between the two models. 
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4.6. Relative intensities of ro-vibrational lines 



We now examine how optical depth effects affect the relative intensities of the CO v=l-0 
ro-vibrational lines. 

As already mentioned, for beams passing close to the nucleus, the ratio of the P2 to RO 
line intensities was measured to be lower than the statistical value of 2 expected under optically 

[)). The East- West averaged P2/R0 intensity ratio was found 



thin conditions (DiS anti et al. 



2001 



to approach the statistical ratio at distances p > 2300 km (> 1.5") in January, and p > 5800 
km (> 5.5") in March-April. This trend can only be explained by opacity effects. As shown in 
Fig. [51 at the temperatures encountered in the Hale-Bopp coma, the opacity of the R0 line is 
smaller than the P2 line opacity (by a factor of 1.6 for January- April) . Figure [131 plots the P2/R0 
intensity ratio derived from our radiative transfer models for the two periods. The on-nucleus 
value of 1.55 found for early March a grees with (although somewhat higher) the measured value 



of ~1.3± 0.1 derived from Fig. Al of 



DiSanti et al 



(|200ll ). not considering the possible slit offset 



in Dec discussed in the previous section (in which case we may conclude that the model slightly 
underestimate optical depth effects). The computed evolution of the P2/R0 intensity ratio along 
the slit agrees with the measurements. As for January, radiative transfer calculations in the 
assumption of an isotropic coma lead to a P2/R0 ratio lower than the observed value (~ 1.3±0.2). 
Our modelling does not predict a significant decrease of the P3/R1 intensity ratio with respect the 
optically thin statistical value of 1.5, because the two lines have approximately the same opacity 
(Fig. [8]). A P3/R1 i ntensity ratio signifi cantly lower than 1.5 was observed, however, but only in 



January (Fig. Al of 



DiSanti et al 



2001 



Overall, there are some indications that our modelling 
underestimates opacity effects, which strengthens our interpretation of the CO spatial profiles as 
being strongly affected by opacity effects. 

In optical thin conditions, the rotational temperature T^ ot describing the population 
distribution in the vibrational excited state can be retrieved from the relative intensities of the 
ro-vibrational IR lines. This temperature is close to the rotational temperature in the ground 
vibrational state T ro t, and probes the gas kinetic temperature when the field of view encompasses 
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essentially the c ollision-domina t ed co ma. Applying the method of rotation diagrams to the CO 
Hale-Bopp data, 



DiSanti et al. 



(J200l|) retrieved T^ ot for v=l, as a function of projected distance 
to nucleus. However, this method is in principle not appropriate since the CO lines are optically 
thick. We applied the same method to synthetic line fluxes, giving the same weight to all lines. 
We found that : i) T^ ot (meas.) < T^ ot , where T^ ot (meas.) is the value retrieved from the rotation 
diag ram; ii) T^ ^{jn€.(is?) and X^ j differ by an amount which depends on the lines used in the 
rotation diagram; iii) T^ ot is larger than T ro t in the ground vibrational state by typically 8-10%. 
For the December simulations performed with T = 60 K (R2-R3, P2-P3 lines), we found for 
the inner l"xl" box: T^ ot (meas.) = 60 K, T^ ot = 66 K; for the January simulations with T 
= 80 K (R0 to R6, PI to P4 lines): T^ ot (meas.) = 65 K, T r v ot = 86 K. The rotation diagram 
method underestimates more significantly T" ot when the R0, Rl and PI to P3 lines are considered, 
because these lines, which are pumped from the lowest energy v = rotational states, are less 



optically thick than the other lines (Fig. [8j). 



DiSanti et al 



(|200ll ) points out that the measured 



temperatures on 24 February 1997 are higher compared with those on most other dates. We 
suggest that this is because the PI to P3 and R0-R1 lines were not used (they were not observed 
on 24 February). In summary, T^ ot values derived from set of lines that include R0-R1 and PI 
to P3 lines likely underestimate (by typically 10-15 K) the gas kinetic temperature T in the 
innermost lines of sight. 



5. Summary 



The high productivity of comet Hale-Bopp in carbon monoxide made possible the investigation 
of the spatial distribution of this molecule in the inner (r = 1000-10 000 km) coma by means 
of spectroscopy from ground-based telescop es. Because there was observational ev idence of an 



extended production of CO in infrared data (jDiSanti et al 



2001 



Brooke et al. 



20031 ) , we performed 



a detailed analysis of the CO data acquired with the IRAM Plateau de Bure interferometer close 
to perihelion (March 1997) for a confirmation of this extended production. This study benefited 
from the simultaneous observations of two CO rotational lines, the combination of both high and 
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low angular resolutions probing coma radii from ~800 to ~20 000 km, and last but not least very 
high signal-to-noise ratios. The results obtained from this study urged us to analyze the CO 
infrared data using radiative transfer codes similar to those used for interpreting the millimeter 
data. We reached the following conclusions: 

1. Both J(l-O) and J (2-1) radio lines exclude pure CO production from an extended source, 
unless the characteristic parent scale length L p of the distributed source was smaller than 
500 km; 

2. Composite models considering half of the CO production from a distributed source cannot 
explain the radio observations, unless L p was < 1500 km. A good fit to the data is obtained 
when CO is assumed to be solely released from the nucleus; 

3. The J (1-0) line data are best explained when considering the radial evolution of the 
temperature in the coma due to photolytic heating. Secure results are obtained for the 

J (2-1) line thanks to the development of a radiative transfer code which accounts for optical 
depth effects in the emergent emission; 

4. The presence of a CO spiralling jet in the Hale-Bopp's coma does not affect much our 
analysis because the millimetric observations spanned a large fraction of the nucleus rotation 
and probed almost the whole coma at projected distances p < 20 000 km; 

5. The observed CO ro-vibrational lines are optically thick. The broad extent of the CO 
brightness distribution in the infrared is mainly due to optical depth effects entering in the 
emitted radiation. Surprisingly, our conclusion is not consistent with the one obtained by 



Brooke et al 



( 2003 ) who performed similar radiative transfer calculations. Comparison of 



model results shows that they quantitatively differ. 

6. The slit-oriented spatial profiles of the CO infrared lines are strongly affected by the CO 
jet, making their interpretation difficult. They are more extended than expected for a pure 
nuclear origin possibly also because the slit was centred on the peak of the dust continuum 
emission which departed from the nucleus position. 



-31 - 



7. The comparison between observed and computed P2/R0 line intensity ratios suggests that 
our modelling underestimates optical depth effects, which strengthens our interpretation of 
the broad extent of the IR emission. 

8. When optical depth effects are present, the rotation diagram method introduces errors in 
rotational temperature determinations. 

In conclusion, our analysis shows that there is no compelling evidence of extended production 
of CO in Hale-Bopp's coma. We anticipate that a similar conclusion should be reached for comet 
Hyakutake from a reanalysis of the observed IR lines. 
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APPENDIX : Visibilities for long-lived parent molecules 

In this Appendix, we demonstrate that, for a parent molecule distribution, the amplitude of 
the visibility V(a) varies with the Mf-radius a according to a -1 , in first approximation. 

Let us start from the general equation of the complex visibility V(u,v) that can be 
approximated to: 

c r+co r+oo 

V(u,v) = - / A(x,y)F(x,y) 

" J -co J — oo 

. 2mu ., , , „. 

x exp( lux + vy)) axay (1) 

c 

in units of line area (Jy km s _1 ). u and v are the coordinates of the baseline vector for two 
antennas in the uf-plane. F(x, y) is the brightness distribution in the plane of the sky. A(x, y) is 
the antenna primary beam pattern of half-power beam width 6>b, approximated by a 2-D Gaussian 
centred on (0, 0) given by: 



a, ^ / 41n(2)(x 2 +y 2 X 
A(x,y)=eM ^ — )■ 

Under optically thin conditions and assuming that the population of the upper level of the 
transition p u does not vary with distance to nucleus, F(x, y) is proportional to f*°° n(r)dz, where 
n(r) is the local density of the considered molecule at distance r from the nucleus. If we assume a 
Haser parent molecule distribution and neglect the photodestruction of the molecule, then: 

1 K 

F(x,y) = K =—, 
\J x 2 + y 2 P 

p is in radians. A is a constant factor equal to: 



-33 - 



_ huA u ip u Q 
167rw cxp A ' 

Q is the gas production rate and v exp is the gas expansion velocity. A is the geocentric 
distance. In polar coordinates in the Fourier (u = a cos a, v = a sin a) and image (x = pcosO, 
y = psin9) planes, Eq.([T]) reduces to: 



Via, a) = / / A(p) exp( pa cos(a — 9)) dOdp. 

v Jo Jo c 



It can be demonstrated that: 



I exp( pa cos(a — a)) dO = 2irJo(- 

Jo c 



c 



(2) 



where Jq is the Bessel function of first kind and zero order. Therefore, Eq.([2]) reduces to: 



. . 2ttcK f+°° A . . T ,2-Kvpa.j 

V{a,a) = / A(p)J ( —)dp. 

v Jo c 

V(a, a) is a real function (a general property of the Fourier Transform of a symmetric function) of 
amplitude V(a) = V(a,a), which can be also written: 



c 2 K 1 f +oc co 
V{a) = °— x - / A{-^)Up)dp. 
a Jq Zirva 

Let us call D the diameter of the antennas. For a/D ^> 0.2 (typically a > 20 m for 15-m 
dishes), the A term in the integral can be considered as unity. As f + °° Jo{p)dp = 1, then: 
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V(a) = hc2A ^Q x I. (3) 

This formula remains valid when photodissociation occurs, provided that a 2> cA/vL, where 
L is the photodissociation scale length. 

In the same units (Jy km s _1 ), the line area for autocorrelation (ON-OFF) spectra which 
corresponds to V(0) is: 



m = (4) 



where N is the number of molecules sampled by the primary beam, which can be computed by 
volume integration. 
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Table 1: Observed intensities of CO lines in single-dish (SD) and interferometric (Int.) modes. 



Line 


Mode 


Beam 


Intensity" 








(Jy km s _1 ) 


J(l-O) 


SD 


41.8" 


10.8 ± 0.43 




Int. 


2.57"x3.58" 


0.93 ± 0.04 






Rl-0 — FSV/ -Pint 


11.6 ± 1.8 6 


J(2-l) 


SD 


20.9" 


82.8 ±0.59 




Int. 


2.0" x 1.38" 


4.65 ±0.15 






= -^Sd/ -Pint 


17.8 ± 3.8 b 



Bockelee-Morvan et al. 



(}2Q09j ). 



a From 

6 Include uncertainties in fl ux calibration: 10% for CO J(l-0), 15% for CO J(2-l) 



( jBockelee-Morvan et al. 



20091). 



Table 2: Modelled PdBI intensities of the CO lines in comet Hale-Bopp and slope of the visibility curves. 



Model 


vco :( *co 


rn 
1 


^exp 


T 

Lip 


-^Int 


-Rmod 


Slope b 




(%:%) 


(K) 


(km s i ) 


(km) 


f \\r k m 
i j y .Eviii 


b ) 




















J(l-0) 


J (2-1) 


J(l-0) 


J (2-1) 


J(l-0) 


J (2-1) 


Nuclear+Thin c 


100:0 


120 


1.05 





0.74 


5.86 


15.0 


14.7 


-1.02 


-1.04 


Nuclear 


100:0 


120 


1.05 





0.70 


4.25 


15.7 


19.1 


-1.07 


-1.20 


Nuclear 


100:0 


rp d 
-L var 


1.05 





1.00 


4.64 


11.6 


18.7 


-0.96 


-1.42 


Nuclear 


100:0 


120 


"^var 





0.78 


4.77 


13.6 


17.2 


-1.04 


-1.21 


Nuclear 


100:0 


rp d 
± var 


V e 
'-'var 





1.10 


4.94 


10.3 


17.7 


-0.95 


-1.50 


Nuclear 


100:0 


T f 

1 var{alt) 


v e 

v var 





1.09 


5.38 


10.7 


16.3 


-0.98 


-1.35 


Extended 


0:100 


120 


1.05 


2000 


0.32 


1.31 


32.2 


55.8 


-1.75 


-1.96 


Extended 


0:100 


rp d 
± var 


^var 


2000 


0.53 


2.27 


19.6 


34.3 


-1.44 


-1.82 


Composite 


50:50 


120 


1.05 


2000 


0.51 


2.92 


20.8 


26.5 


-1.26 


-1.35 


Composite 


50:50 


120 


1.05 


5000 


0.44 


2.60 


23.3 


27.5 


-1.22 


-1.26 


Composite 


50:50 


rp d 
1 var 


v e 

'-'var 


2000 


0.84 


3.89 


13.0 


21.3 


-1.07 


-1.51 


Composite 


50:50 


rp d 

- 1 var 


v var 


5000 


0.72 


3.48 


14.3 


21.9 


-1.03 


-1.40 


Composite 


10:90 


rp d 

± var 


v e 

v var 


5000 


0.37 


1.76 


25.4 


37.7 


-1.33 


-1.53 



-41 - 



Footnotes to Table 2 : 

a Relative contribution (in percent) of the nuclear CO production rate Qqq and extended CO 
production rate Qqq *° ^ ne total CO production rate Qco = Qco + Qco = ^-1 xlO 30 s _1 . 

b Slope (i.e., power index) of the variation of the visibility amplitude with uf-radius between 20 
and 150 m. 

c Opacity effects are not considered in radiative transfer calculations. 

d Variable kinetic temperature T var : T = 20 + 50 x log(r[km]/100) K for 100 < r < 10000 km, 
T = 120 K for r > 10000 km, T = 20 K for r < 100 km. 

e Variable expansion velocity v var : v cxp = 0.9 + 0.15 x log(r[km]/1000) km s" 1 for r > 1000 km, 
Vcxp = 0.9 km s" 1 for r < 1000 km. 

f Variable kinetic temperature T var(alt) : T = 20 + 50 x log (r [km]/ 100) K for 1500 < r < 10000 km, 
T = 120 K for r > 10000 km, T = 79 K for r < 1500 km. 
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Table 3: CO IR line fluxes in March 1997. 





T 1 Tl 
J-J111C 


I mm in f^irt" 

V^wllllllCLlu 


vco 


U 1 J L L 


p. b,c 
± inner 


J? b,d 
r outer 


F * IF- 

1 outer 1 1 inner 








(s" 1 ) 


(") 


(W m~ 2 ) 


(W m" 2 ) 




3.80 March UT 


R2 


Measured b 






43 ± 1.5 


46.8 ± 2.9 


1.09 ± 0.08 






Panel 3/Iso b 


2xl0 30 





43.1/45.5 


37.8/36.6 


0.88/0.80 






Panel 3/Iso b 


2xl0 30 


-1.2 


40.7/39.1 


37.4/35.8 


0.92/0.92 






Panel 3/Iso b 


2xl0 30 


-2.4 


31.5/28.6 


33.6/33.7 


1.07/1.18 


5.99 March UT 


R6 


Measured^ 






20.6 ± 1.6 


21.4 ± 1.4 


1.04 ± 0.11 






Panel 10/Iso 6 


2xl0 30 





36.9/34.9 


31.2/26.1 


0.85/0.75 






Panel 10/Iso 6 


2xl0 30 


-1.2 


26.8/29.2 


29.8/25.6 


1.11/0.88 






Panel 10/Iso fe 


2xl0 30 


-2.4 


20.9/20.6 


26.7/24.1 


1.28/1.17 



a Slit offset in declination. The slit is oriente c 
b Measured intensities are from Brooke et al. 



last- West. 



(120031 ). Modelled intensities were computed 



using the brightness distributions of Panels 3 and 10 of FigdU for 3 and 5 March, respectively 
(first number in columns 6-8) and an isotropic model (referred to a s Iso, second number 



2003] ), for 3 and 5 March, 



in columns 6-8). A seeing of 1.78" and 1.7" (IBrooke et al. 
respectively. 

c Line intensity in a 2"xl" box in RAxDec centred on the peak of CO brightness along the 
slit. 

d Line flux in East- West averaged 5"xl" boxes beginning 3 " from peak. 
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Fig. 1. — ON-OFF and interferometric spectra of the H?S 2^-2^ line at 2 16.7 GHz obtained 



20071 ). Synthetic spectra 



in comet Hale-Bopp on 13 March 1997 with the PdBI (IBoissier et al. 
computed using the v var law for the expansion velocity are shown with dotted lines. The 
angular resolution is 22.2" and 1.7", for the ON-OFF and interferometric beams, respectively. 
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Fig. 2. — x 2 — (Robs — -RMod) 2 /^-Robs as a function of the CO parent scale length, where 
.Robs and -RMod are the observed and modelled intensity ratios R— Fsr>/Fi nt , respectively, 
and Ai?obs is the error bar on -Robs- Blue (respectively red) symbols are for the CO J(l- 
0) (respectively J(2-l)) lines. Results plotted in the right part of the plots (positive scale 
lengths) are for pure extended distributions (i.e., Qco'-Qco = (0:100)). Results for composite 
distributions {Qco'-Qco = (50:50)) are shown in the left part of the plot (negative scale 
lengths). The grey region corresponds to < % 2 < 1. The four quadrants correspond to 
different gas temperature and velocity laws indicated in the bottom right corner. 
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Fig. 3. — Amplitude of the visibility as a function of wf-radius for the CO J(l-O) line 
(top) and the CO J (2-1) line (bottom). Measurements are shown by dots with error bars 
which correspond to statistical noise uncertainties, except for the ON-OFF measurements 
where calibration uncertainties (10% and 15% for J (1-0) and J (2-1), respectively) are added 
quadratically. The different curves correspond to different models: nuclear production with 
variable gas temperature in the coma (plain line), nuclear production with constant temper- 
ature (dotted line), combined nuclear /extended productions in ratio 50:50 with L p = 2000 
km and variable temperature (dashed line), combined nuclear/extended productions in ratio 
50:50 with L p = 5000 km and variable temperature (dotted-dashed line). All models consider 
isotropic outgassing at a velocity described by t> var and a total CO production rate Qco = 
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Fig. 4. — Log-Log plot of the amplitude of the visibility as a function of wf-radius for the 
CO J (1-0) line (top) and the CO J (2-1) line (bottom). Model parameters are the same as 
for Fig. [31 
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Fig. 5. — Amplitude of the visibility as a function of wf-radius. CO J (1-0) (respectively 
J (2-1)) measurements are shown with open (respectively filled) symbols with error bars. 
The different curves correspond to calculations for a CO nuclear production and different 
temperature laws. Other model parameters are the same as for Fig. [3J 
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Fig. 6. — Same as Fig. [2} using the T var ^ a i t ) and v var laws 
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Fig. 7. — Q-curves generated with the isotropic model of CO v = 1-0 line emission. CO is 
assumed to be released from the nucleus. Lines used are R2, R3, P2 and P3 (Dec. data), 
RO, Rl, PI, R2, P2, R3, P3, R4, R5, R6 (Jan. data). The aperture is set to l"x 1". These 



Q-curves can be directly compared to Fig. 7 of iDiSanti et al.l (120031 ). Plain symbols: results 
with all processes included (optical depth effects, full rotational and vibrational excitation 
model). Open symbols: reference Q-curve, i.e., model results under optically thin conditions, 
assuming Boltzmann rotational level populations and ineffective vibrational collisional relax- 
ation. Dotted curve: optically thin Q-curve scaled to the value of the optically thick Q-curve 
at the nucleus position. Calculations were performed with T = 60, 80 K, and v exp = 0.9, 1.0 
km s _1 , for December 11, 1996, and January 21, 1997 conditions, respectively 
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Fig. 8. — Mean opacity of observed CO IR lines for conditions of early-March 1997 (A = 
1.47 AU, Qco = 2 x 10 30 s" 1 , T = 90 K, v exp = 1.1 km s _1 ). Calculations are for a circular 
aperture of 1" diameter centred on the nucleus. The mean opacities of the 115 and 230 GHz 
lines in the same FOV are plotted, for comparison. 
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Fig. 9. — Same as Fig. [7] with gas acceleration considered in radiation transfer calculations. 
Apparent CO production rates were computed assuming a constant outflow velocity v exp = 
0.9 and 1.0 km s _1 , for 11 December 1996, and 21 January 1997, respectively. 
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Fig. 10. — Same as Fig.[7Jfor January 21, 1997, considering both pure nuclear production and 
combined (50:50) nuclear /extended distributions: nuclear (filled dots), combined with L p = 
1500 km (stars), 2500 km (empty triangles) and 5000 km (crosses). The reference Q-curve 
is shown by empty circles connected by dashed lines, as in Fig. [7J Scaled reference Q-curves 
are shown with dotted lines. Other model parameters are the same as for the calculations 
shown in Fig. [7J 
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Fig. 11. — Time evolution of the brightness distribution (not corrected for seeing) of the 
v = 1-0 R6 CO line f o r geo metrical conditions of 1-5 March 1997. The coma model of 
ockelee-Morvan et all (120091 ) is used (jet model (3); Q co = 2 x 10 30 s" 1 ). The gas coma 
temperature is taken equal to 90 K. Panels 1 to 12 cover one rotation period P = 11.35 h 
with a time step of P/12. Line fluxes were computed for 0.2x0.2" pixel size. Contour levels 
are in logarithmic scale and spaced by 0.15. Maximum fluxes on the maps are between 9.2 
and 9.6 10 -18 W m -2 . The nucleus North pole (N) and the Sun (S) directions are indicated. 



Celestial East is on the left. 
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Fig. 12. — CO v = 1-0 R6 spatial profiles derived from the brightness maps (Panels 3, 4, 
10, and 11) shown in Fig. [TT] convolved with the PSF. Extractions are parallel to the RA 
axis and centred at SDec = 0" (plain line), <5Dec = -1.2" (dashed-dotted line) and 5Dec = 
-2.5" (dotted line). Intensity has been integrated over rectangular boxes of 1" along Dec 
axis and 0.2" along RA axis. The dashed line shows a column density profile convolved 
with the PSF. The seeing was taken equal to 1.7 ". Peak intensities (units of 10 -17 W m -2 ) 
are 3.5, 3.2, 2.4 (Panel 3), 3.5, 3.4, 2.8 (Panel 4), 3.4, 2.7, 2.1 (Panel 10), and 3.3, 2.6, 1.9 
(Panel 11) for <5Dec = 0, -1.2", and -2.4", respectively. 
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Fig. 13. — Ratio of the P2 to RO line intensities along a slit aligned along the RA axis and 
centred on the nucleus. Results from radiative transfer calculations in an isotropic coma are 
shown in dashed (parameters of 21 Jan. 1997) and plain (early March 1997) lines. Results 
obtained with the jet model (phase of panel 10) are shown in dotted lines. The ratio was 
computed using PSF-corrected intensities over rectangular boxes of 1" along Dec axis and 
0.2 " along RA axis. 



